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ABSTRACT 



The sensitivity of the XR-3 roll behavior in calm 
sea on sidewall parameters and force and moment 
calculations is investigated with the Loads and 
Motions computer program and compared with 
experimentally measured data. Propulsion and rudder 
subroutines and added mass computation are reviewed 
and modified. Recommendations for improved simulation 
of the XR-3 roll behavior are given. 
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I. introduction 



A. Background 

In the past several years the Surface Effect Ship (SES) 
has received increased attention in the United States Navy 
and detailed investigations of its sea performance have been 
carried out. The two major categories of SES are the Air 
Cushion Vehicle (ACV) and the Captured Air Bubble Ship 
(CAB) . The CAB-type of craft has become of primary interest 
to the United States Navy and two of the constructed 
testcrafts are the Bell Aerospace Systems 100-B and the XR-3 
of approximately 100 and 3 tons of displacement, 
respectively. In 1976 the U.S. Navy Surface Effect Ship SES 
100-B of the CAB-type established a world speed record for 
surface type ships of 89.4 knots. In April 1976 the SES 
100-B launched a missile at 60 knots. The launch was 
successful and the missile hit its target five miles away. 



In order to investigate the characteristics of the 
CAB-type SES under any design and operating conditions 
without the costly need of modifying the actual craft, a 
Loads and Notions digital computer program was developed by 
Oceanics, Inc. [Ref. 1]. In 1970 the XR-3 testcraft was 
delivered to the Naval Postgraduate School and Leo and 
Boncal converted the 100-B Loads and Motions program to 
represent the XR-3. Since there were substantial design 
differences between these two ships programming 
modifications were required in certain subroutines. 
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B 



Objectives 



The purpose of this thesis is to inve 
of safe maneuverability of these high s 
their roll behavior in a turn maneuver 
objective the sensitivity of the simu 
behavior in turn maneuvers as it is repres 
and Motions Program is investigated for i 
changes in sidewall characteristics, 
modifications in the added mass compu 
parameters and force calculations as well 
rudder input parameters are investigated, 
the computed performance is obtained 
experimentally measured steady state roll 



stigate one aspect 



peed 


ships. 


i. e. 


. To 


fulfill 


this 


lated 


XR-3 


roll 


ented 


by the 


Loads 


ts dependence on 



The effect of 
tation, sidewall 
as propulsion and 
An evaluation of 
by comparison with 
conditions. 
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II 



INIT IAL RE.^ARKS 



A. PROBLEM OF REPEATABILITY 



It has been this author's experience as well as all 
previous investigators at the Naval Postgraduate School 
using the XR-3 Loads and Motions Program that there have 
been no observed problems related to repeating the results 
for a given input condition. However, during the author's 
first studies of the Loads and Motions Program as listed in 
Menzel's thesis [Ref. 2] an attempt was made to duplicate 
some of the simulation runs given in that reference work. It 
was found necessary as a first step to restore the 

simulation program on the IBM 360/67 computer from a card 
source onto a disk. ^ Using the input conditions stated in 
Ref. 2 it was impossible to obtain identical time 

histories, e.g. for turns. Each time subroutine SIDEWALL was 
read into the computer together with the data input deck, 
thereby overriding the stored program, an error message in 
the printed output indicated that a division by zero 
occurred in this subroutine. The error message was not 
generated when the stored program which contained an 
identical version of subroutine SIDEWALL was called. A 
closer look revealed that in the SIDEH ALL- version given in 
Ref. 2 the variable PBAR (plenum pressure) was neither 



1 It should also be noted that results obtained in Ref. 2 
used the IBM 360/67 Release 20.6 system while the results 
obtained in this work used rhe 21. 8B Release with HASP 
installed. 



defined nor transferred by a COMMON-statement and therefore 
a default value of zero was used during the computation of 
PBHEAD (see Appendix C, SDWL 0440). The missing statement 

PBAR = PB - PINF 

was inserted to subroutine SIDEWALL. With this correction 
the new results were slightly closer to those given in 
Ref. 2. 



The table given below compares output values for roll 
and pitch angle in a turn at 20 kn with constant thrust and 
a 35 degrees rudder step input under calm sea conditions. 



20 k n tur n a t 35 deg ru d der angl e 



without with 
PBAR-card PBAR-card Ref. 2 



relative 
de via tion 



Roll angle (deg) 
first peak 
avrg at 20 sec 

Pitch angle (deg) 
first peak 
avrg at 20 sec 



7.44 7.39 

2.65 2.77 



1.34 1.33 

0.99 0.94 



7.2 2.6 % 

3.1 1 1 . 9 7o 



1.2 9.3 % 

0.8 1 4 . 9 -;5 



Although the differences between columns 2 and 3 in TABLE I 
might be considered to be small as far as magnitudes are 
concerned , the relative deviation is up to 15 % as shown in 
column 4. But what is more important is the fact that the 
roll and pitch angle responses for the uncorrected (plots 1 
and 2 in Appendix A) as well as for the corrected program 
(plots 3 and' 4) show an unstable craft behavior for t>40 sec 
where the pitch angle increases rapidly while the roll angle 
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decreases rapidly approaching the zero degree value. This 
unstable condition is probably due to the 35 degrees rudder 
step input at 20 kn, as it is used in Menzel's study, which 
is physically unreal since it could never be introduced to 
the XR-3 testcraft in an actual run at that speed. 
Therefore rudder angles of up to 15 degrees introduced at a 
rate of 5 deg/sec will be used throughout this work, as it 
is also done on the testcraft. 



also, regarding plots 1-4 again and 
whether they really reflect the actual cr 
should be noted that a short time history 
Ref. 2) could possibly result in wrong con 
the unstable condition is not evident at t 
time. Therefore, time histories of up to 
will be shown throughout this study. 



not consi 
aft behavi 
( t<20 s 
elusions b 
hat insta 
40 or 50 s 



dering 
or , it 
ec in 
ecause 
nt of 
econds 
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B. STEADY STATE CONDITIONS 



Due to the 


change in 


subroutine SIDEWALL 


as mentioned 


the preceding 


section 


the steady 


state 


conditions 


straight runs 


in calm 


water have 


been reestablished 


various speeds 


and are listed in Table 


II. 








Table II 








Steady State Conditions 




Speed 


Draft 


Pitch 


Plenum 


Thr us t 






angle 


pressur 




(kn) 


(in) 


(deg) 


(psf) 


(lb) 


10.0 


8. 17 


1.62 


23 . 93 


400 .61 


12.5 


7.03 


1.11 


24 . 84 


349.42 


15.0 


6.74 


0.84 


24 . 84 


335.71 


17.5 


6.41 


0.63 


24 . 34 


346 .45 


20.0 


6.12 


0.48 


24 . 84 


373 .44 


22.5 


5.87 


0.36 


24.84 


411.53 


25.0 


5.66 


0.29 


24 . 84 


458.62 


27.5 


5.48 


0 . 25 


24 .84 


513 .31 


30.0 


5.34 


0.26 


24.84 


574 .43 



These steady state conditions have been established by 
first executing the XR-3 Loads and Motions Program with 
constant speed for 40 seconds and then repeating the run 
keeping the evaluated thrust constant. A comparison of 
these steady state values with those given in Ref. 2 shows 
differences especially in the lower and upper speed range of 
again up to 15 % the cause of which could not be suspected 
anywhere else since the simulation program has not been 
changed since last used by Menzel. 
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c. 



SIGN CONVENTIONS 



The sign conventions used in the XR-3 Loads and Motions 
Program are identical to those used in Report 71-84 by 
Oceanics Incorporated [Ref. 1]. A right handed coordinate 
system is applied to the craft with the positive X-axis 
being measured forward and the lateral Y-axis being measured 
positive to starboard. The vertical Z-axis is measured 
positive downward. Identical signs are used for respective 
forces along those axes. 

The signs of the angles are defined in the following 
manner : 

pitch angle being positive upward 

(boat noses up) 

roll angle being positive to starboard 

(boat rolls to starboard) 

yaw angle being positive to starboard 

(boat turns to starboard) 

rudder angle being positive to port 

(right rudder, boat turns to port) . 

Zero pitch and roll angle are referenced to the X-Y 
plane parallel to the water surface. 
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III. INVESTIGATION OF SIDEWALL EFFECTS 



A. ADDED MASS EFFECT 



1 . Theory 



A basic element in representing the hydrodynamic 
forces by use of slender body theory are the two-dimensional 
sectional values of lateral added mass A22 of the sidewalls, 
which are necessary for the determination of the lateral 
forces, as well as the two-dimensional sectional vertical 
added mass A33 which is also used in determining the roll 
moment. These two-dimensional added mass values are given in 
the 'Surface Effect Ships Aero/Hydrodynamics Technology 
Design Manual' [Ref. 3 ] as 



and 



A22 



C *RH0*PI*D«D/2 . 0 



A33 = C *RHO*B*3/8.0 

V 



where C is the lateral added mass coefficient 
h 

C is the vertical added mass coefficient 

V 

RHO is the specific density of sea water 
D is the local draft 

and B is the width at the local sidewall waterline. 
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were 



The C and C values originally selected in Ref. 1 

h V 

those corresponding to high speed with C =0.4 and C =1.0. 

h V 

But in order to account for variations in sidewall shape, 

and their influence on the effective lateral added mass, the 

results of the research work described in Ref. 3 led to an 

average value of C =0.8. This new coefficient improved the 

h 

agreement between theoretical and experimental data as 

stated in Ref. 3 . In the present XR-3 Loads and Motions 

Program a value of C =0.4 is used in accordance with Ref. 1 

h 

and the effect of letting C =0.8 will be investigated in 

h 

later parts of this thesis. 



Considering the equation for the vertical added mass 
A33, a major difference has been found between that one 
given in Ref. 3 and the one used in the XR-3 Loads and 
Motions Program (e.g. Ref. 2), where in accordance with 
Ref. 1 (1971) there was an additional Pl-factor in the A33s 
equation. The letter • s' indicates that the computation of 
the added mass is done at the craft's stern. Since Ref. 3 
(1976) reflects the experimental and theoretical work being 
done based on Ref. 1, the A33s equations in subroutines 
SIDEWALL and SIDETAB have now been changed to the new form 
given above. 
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2 . 



Ne w . Ste ady State Conditi ons 



The change in the simulated XR-3 performance due to 



the reformulated 


A33s computation can be 


observed 


from 


TABLE III giving the new steady state 


conditions 


for 


straight runs in 


calm sea 


for various speeds which 


are 


identical for 


both (0. 


4 and 0.8) lateral added 


mass 


coefficients . 














T^LE_III 








Steady 


state conditions 










(new A33s) 






Speed 


Draft 


Pitch Plenum 


Thr usr 








angle pressure 




(kn) 


(in) 


(deg) (psf) 


(lb) 




15.0 


6. 80 


0.86 24.84 


336.30 




17.5 


6.46 


0.65 24.84 


347.23 




20.0 


6. 17 


0.49 24.84 


374. 38 




22.5 


5.92 


0.37 24.84 


412. 62 




25.0 


5.71 


0.30 24.84 


459.90 




27.5 


5.53 


0.27 24.84 


514. 89 




30.0 


5.40 


0.28 24.84 


576.66 




Comparing these 


steady 


state values with 


those previously 


given in TABLE II in Sect. 


II. 3 one finds in 


general for 


all 


speeds that 










- draft 


has increased by about 0.05 


in (0.8 %) 




- pitch 


angle has 


increased by 0.02 


deg (2 to 7 


%) 


- thrust 


has slightly increased by less than 0.5 


or 

/O • 
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3. 



Effect 1 r a T urn m an euve r 

The effect of the change in the formulation of the 
vertical added mass in a simulated turn maneuver has been 
studied next. After a 5 sec straight run at 20 kn in calm 
sea a port rudder deflection was introduced at a rate of 
5 deg/sec and then kept fixed at 15 degrees resulting in a 
turn to port. The deadrise forces are computed at the 
transom. In TABLE 17 are shown the steady state conditions 
and the roll moments contributed from the bow seal (FKBS) , 
stern seal (FKSS ) , sidewalls (FKSW) , rudder (FKRUD) , 
propulsion system (FKP) , aerodynamics (FKAED) and plenum 
pressure (ABPB*PHI*Z) . 



T ABL E IV 



Stead 7 


state conditions in 20 kn 


turn 




1 5 deg port 


rudder angle 








old A3 3s 


new A33s 


dev iati 


on 


pitch angle (deg) 


0.48 


0.52 


8.3 


% 


roll angle (deg) 


2.0 


1 .98 


- 1.0 


% 


draft (in) 


5.36 


5.94 


1.4 


% 


speed (kn) 


19.25 


19 .27 


0. 1 


% 


moments: 










FKBS 


-391.1 


-336.3 


- 1.2 


% 


FKSS 


- 3.6 


- 3.6 


0.0 


% 


FKSW 


191.0 


202.0 


5. 8 


% 


FKRud 


- 97.5 


-110.0 


-12.8 


% 


FKP 


0.2 


0.2 


0.0 


% 


FKAed 


- 54.4 


- 53.5 


1.7 


% 


ABPB*PHI*Z 


355.4 


351 .9 


1.0 


% 
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From TABLE IV it can be seen that, due to the 
modified A33s computation, the counteracting roll moments 
from the sidewall (FKSW) and rudder (FKRtID) changed by about 
6 and 13 respectively, resulting in only a little effect 
in generating a tendency toward a smaller roll angle which 
is directed out of the turn. Other major contributions to 
the roll angle are from the bow seal (FKBS directed into the 
turn) and from the plenum pressure (ABP3*PHI*Z) directed out 
of the turn. The corresponding plots (plots 5-8 for old 
A33s, plots 9-12 for new A33s) show that with the new 
equation for A33s the pirch and roll angle response curves 
are less damped than they were using the former equation for 
vertical added mass calculation as would be expected. 

It should be noted that the roll moment due to 
rudder is fairly large and changes quite significantly while 
the moment contributed by the propulsion is very small and 
constant. The cause of this will be investigated in 
chapter IV. 
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V 



B. DSADRISE ANGLE 



1 . B ack g round 

Surface Effect Ships (SES) of the captured air 
bubble type (CAB) have rigid sidewalls immersed into the 
water, thereby - together with the bow and stern seal - 
capturing the air in the planum chamber and reducing air 
leakage. Since the sidewalls of the XR-3 are not uniform in 
cross-sectional shape throughout their length but are curved 
on the outboard side near the bow similar to a single hull 
ship, the expression '•sidehulls” would be more appropriate, 
as can be found in modern literature. 

The importance of the correct understanding of the 
effect of the sidewalls on the craft's performance is 
readily seen .from a report on SES Research with the XR-1 
testcraft [Ref. 5]. As reported, the first test series were 
carried out successfully with a 45 degrees sidewall deadrise 
angle, which is the angle between the horizontal plane and 
the ship's outer sidewall surface. In 1964, after a 
modification of the craft to 60 degrees deadrise angle and 
articulating seals, an unstable roll condition occurred 
during a turn at 35-37 knots, a maneuver performed many 
times before. The testcraft heeled out of the turn , nosed 
down due to the retraction of the bow seal, continued its 
outward roll motion and then flipped over. The report, 
however, does not mention the rudder action generating the 
turn. After this accident, the deadrise angle had been 
changed back to 45 degrees. From this experience the 
surface effect ship's roll stability should be expected to 
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be sensitive to the deadrise angle and an investigation of 
this sensitivity is reported in the following sections. 

2 • Force s d ue to Deadri se Angle 

The deadrise angle, as pointed out by Ref. 3 and 4, 
has a major effect on the craft's roll motion and therefore 
is an important design consideration. During a turn the 
craft is desired to heel inboard, thereby minimizing the 
apparent transverse acceleration (a coordinated turn). For 
typical 3ES designs, it is not possible to achieve perfect 
coordination, but it is possible to achieve nearly flat 
turns if the craft is designed properly. The principal 
factors affecting roll stability at speed are 

* sidewall geometry 

* seal characteristics 

* vertical location of the center of gravity. 

The forces acting at the center of gravity station of the 
starboard sidewall in a port turn are shown in Figure 1. 
The centrifugal force generated acts away from the center of 
the turn. This force must be counteracted by hydrodynamic 
forces which are created by yawing the ship into the turn. 
To an observer aboard the craft, it would appear as if the 
craft is sideslipping. Wave buildup on the outboard sides of 
the sidewalls increases the pressure there, while only small 
pressure changes occur on the inboard side of the opposite 
sidewall generating a small (negligible) vertical force. 
The principal force arises from the outboard side of the 
sidewall. The direction of this force depends on the slope 
of the deadrise angle. If the resultant force passes above 
the center of gravity (solid line in Fig. 1) a restoring 
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moment results, tending to heel the craft into the turn 
thereby improving its roll stability. If the deadrise angle 
is chosen to be larger (e.g. 70 deg, shown dashed), the 
resultant force passes below the center of gravity, 
resulting in a moment tending to heel the craft out of the 
turn . 



Thus, once the SES geometry and the exact vertical 
location of the center of gravity are known, an approximate 
check for stability in a turn can be made. 




Figure 1 
(C. F. 



- GENERAL ROLL STABILITY CHECK IN A PORT TURN 
CENTRIFUGAL FORCE, C. G. CENTER OF GRAVITY) 
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XR-3 Sidewall Geome try 



From the construction data entered into the sidewall 
table (subroutine SIDETAB) Figure 2 has been drawn showing 
the cross sections at several selected stations of the 
starboard sidewall. Since in most simulation runs (which 
will be discussed later) draft was less than 7 inches, the 
deadrise angles for all stations immersing into the water 
have been calculated by straight line approximation and are 
given in TABLE V. The stations are counted from bow 
(station 0) to stern (station 28) and are not equally spaced 
for station numbers less than 11. The center of gravity has 
been located experimentally at station 14. The average 
deadrise angle of all stations listed is 64.3 degrees. 

TAB LE V 

Deadri se angl es (deq) fo r all st ati ons 

S tra igh t line a ppro ximati on fo r 7 * ' d raft 



ation 


Deadrise angle 


Station 


Deadrise angle 


5 


67.6 


17 


62.1 


6 


70.8 


18 


64 .0 


7 


57.9 


19 


65.3 


8 


54. 4 


20 


66 .7 


9 


53.7 


21 


68. 1 


10 


53.7 


22 


6 . 96 


1 1 


53.7 


23 


70.3 


12 


55.6 


24 


7 1.1 


13 


57. 3 


25 


72.4 


14 


58. 5 


26 


73.8 


15 


59.7 


27 


76 . 9 


16 


60.9 


28 


78.8 
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Figure 2 - XR-3 STARBOARD SIDEWALL CROSS SECTIONS 
(STATIONS AS INDICATED, UNITS IN FT.) 
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4 . Effect o f Deadri se Angle 

From the previous discussion the surface effect 
ship's stability behavior is expected to be sensitive to the 
sidewalls' deadrise angle. In the XR-3 Loads and Motions 
Program listed in Ref. 2 the deadrise angle at the transom 
(78.8 deg) had been used for the computation of the deadrise 
forces. This calculation was implemented into subroutine 
SIDEWALL using the following statements; 

FYH (J)=-A22S*U* (V+XSS*R-ZS*P) 

# 

• 

CTNDR=0.0 

IF (DSS. LE.0.0) GO TO 22 

CTNDR=(ES-3B( 1) ) /DSS 

IF (THETA.lt. 0 .0) CTNDR=0. 39391 
22 CONTINUE 

FZHOLD (J) =FZH (J) 

FZHDRP (J) =PM1*FyH (J) *CT ND R* PROM 0 1 

FZH (J) =FZH (J) +FZHDRP(J) 

• 

FK= (FZH (2)-FZH(1) ) * YSW + FKD-FY+ZS 

The projection of the lateral force ?YH(J) at each sidewall 
(J=1 for port, J=2 for starboard side) produced the 
respective component of vertical force FZH, where FZHDRP is 
the projected force and CTNDR is the cotangent of the 
deadrise^ angle. The roll angle had been already taken into 
account in the computation of DSS (draft at the stern) . The 
factor PfiCMOl simply provided a means to arbitrarily change 
the vertical projected force for studies undertaken in 
Ref. 2. The factor PM1 introduced rhe necessary sign change 
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in the sidewall force computation being dependent on the 
craft's side (FM1=-1 for port, PM1=+1 for starboard side). 
The roll moment FK was partly determined from the vertical 
forces of both sidewalls. 

Comparing the deadrise angles given in TABLE V it is 
seen that the transom deadrise angle of 78.8 degrees is 
unique and the largest along each sidewall and no reason 
could be found why just the transom deadrise angle should be 
used for the vertical force computation. This fact rather 
came up when Leo and Boncal (Ref. 8) scaled down the 
simulation program of the 100-B testcraft to create the XH-3 
Loads and Motions Program. For the 100-B testcraft the 
sidewall cross sections are uniformly shaped throughout most 
of its length. This uniformity does not occur for the XR-3, 
as shown in Figure 2. Therefore it has been supposed that 
another deadrise angle which is more representative for all 
angles existing along the XR-3 sidewalls, e.g. the deadrise 
angle at the center of gravity location (station 14), could 



possibly be 




more appropriat 


e for 


the vertical 


force 


computation . 


To 


investigate the 


effect 


of this 


suppos 


it ion. 


turn maneuvers 


in calm sea at 


20 knots and various 


rudder 


angles have been 


simulated for 


both 


deadrise 


angle 


s and 


lateral added 


mass coefficients. 














TABLE 


VI 









Steady S tate Ro l l Angle (de^) a t 20 icn 

for vari ous Rudder A n gles ( deg) 



Rudder 


Deadrise 


force computation at 


Experimental 


angle 


transo m 


center 


of gravity 


testcra ft 




C =0. 4 


C =0.4 


C =0.8 


data (Ref. 6) 




h 


h 


h 


5 


0.52 


0.32 


0.29 


0.09 


9 


1. 07 


0.78 


0.72 


0.28 


1 2 


1. 50 


1.15 


1.11 


0.58 


15 


1.98 


1 . 53 


unstable 


1 . 36 
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Table VI shows the steady state values for the roll 
angle which are compared with those measured experimentally 
in 1974 and documented in Ref. 6. Since the XR-3 Loads and 
Motions Program under investigation represents the XR-3 
craft configuration as it existed in 1974, Ref. 6 may serve 
to check whether the results produced by certain changes in 
the simulation program are favorable or not. It is not 
anticipated to exactly match the results obtained from the 
simulated runs to the measured values since for the 
testcraft data the following precisions for the measuring 
devices are stated in Ref. 6 ; pitch and roll angle ±0.5°, 
rudder angle ±1.0° and speed ±1 kn. From TABLE VI it is 
obvious that the steady state values for roll angle using 

C =0.4 as the lateral added mass coefficient and the center 
h 

of gravity deadrise angle for the vertical force computation 
are closer by 25 to 40 % to the measured angles than are the 
corresponding roll angles considering the transom deadrise 
angle . 



Osing the center of gravity deadrise angle and the 

lateral added mass coefficient C =0.8 as suggested in 

h 

Ref. 3, the agreement in steady state roll angles could be 



improved by another 3 to 10 ^ for rudder angles up to 

12 degrees. This improvement has also been reported in 

Ref. 3. Comparing the roll angle responses for turns 

generated by a 12 degrees rudder angle to port with C =0.4 

h 

{plots 21 and 22) and C =0.8 (plots 23 and 24) it is found 

h 

that the smaller lateral added mass coefficient generates a 



little higher peak roll angle (1.470) than does the larger 
coefficient (1.27®) . Both responses show the same number of 
oscillations (eleven cycles) until they die out. The larger 
lateral added mass coefficient generates a small negative 
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roll angle shortly after the port rudder motion has been 
introduced. However, for a 15 degrees port rudder angle the 
simulation program showed an unstable response with 
increasing amplitudes of oscillation for pitch and roll 
angle (see plots 17 and 18) . From this the lateral added 
mass coefficient being 0.8 does not seem to be appropriate 
for the XR-3 simulation if both sidewalls are considered for 
the vertical deadrise force computation. The effect of 
using both lateral added mass coefficients will be 
considered again in Sections III.D and E. 



From TABLE 7II it can be seen that the roll moments 
contributed by the bow seal (FKBS) , plenum pressure 
(ABPB*PHI=«'Z) , rudder (FKROD) and sidewalls (FKSW) change the 
most. Comparing the effect due to the change from transom to 
center of gravity deadrise angle it is found that 

- bow seal inward effect decreases 

- plenum pressure outward effect decreases 

- rudder inward effect decreases 

- sidewall outward effect increases 



= outward roll angle decreases. 



The time 
15 degrees port 
transom deadrise 
as plots 13-16 
plots the effect 
gravity deadrise 



histories for simulated turn maneuvers with 
rudder angle using center of gravity and 
angle (TABLE VII) are shown in Appendix A 
and 9-12, respectively. Comparing these 
of the change from transom to center of 
angle and C =0-4 can be observed as 



29 



- pitch angle response being more damped 



- roll angle response being more damped 

- pitch and roll rate responses being more damped 
and having smaller peak values. 



TABLE VII 



Steady State Values for selected 


Deadrise 


A n qles 


20 


kn Turn, calm Sea 


, 15 deq port rudder 


angle 


Deadrise Angle Transom 


C. of G. 


arbitrarily chosen 


(degrees) 78.8 


58.5 


31.6 


28.7 


Pitch angle 


(deg) 0.52 


0.45 


0.52 


0 .64** 


Roll angle 


(deg) 1.98 


1.53 


0.31^ 


- 0.33* 


Draft 


(in) 5.94 


5.81 


6.07 


6.52 


Speed 


(kn) 19.27 


19.37 


19.1 


18.7 


FYSW 


-824.0 


-821.9 


-810.0 


-801 .0 


FYRUD 


40.0 


32.7 


39.5 


55.8 


FYAED 


- 19.5 


- 20.0 


- 19.1 


- 17.6 


R*V*AM 


803.5 


809.2 


789.6 


762.8 


FKES 


-386.3 


-336.4 


- 68.4 


71 .6 


FKSS 


- 3.6 


- 2.7 


- 0.5 


0 .5 


FKSW 


202.0 


212.1 


175.0 


184.0 


FKRUD 


-110.0 


- 89.9 


-108.4 


-153.0 


FKAED 


- 53. 5 


- 56.0 


- 52.4 


- 46.2 


FKP 


0. 19 


0.14 


0.03 


- 0 .03 


ABPB*PHI*Z 


351.9 


273.5 


54.6 


- 57.2 


FYH (P/S) 


-3 9/- 147 


-45/-126 


-78/-95 


-107/-88 


FYD 


-637. 8 


-651.4 


-637.3 


-607 .0 


Note : * = 


still decreasing 








** = 


still increasing 


, not quite 


steady state 
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The reason for the sidewall roll moment and thereby 
the net roll effect not changing more pronounced (TABLE VII) 
is that not only the lateral force FYH must be considered in 
order to determine the force F effecting roll stability 
(Sect. III.B.2, Fig. 1) but also the lateral drag force FYD 
has to be taken into account. Considering the XR-3 geometry 
using the center of gravity deadrise angle, as it is shown 
in solid lines in Fig. 3, the lateral force FYH can be seen 
to generate the vertical projected force FZHDRP and the 
force F, which is directed well above the center of gravity. 
But since FYD has to be added vecrorially to the force F and 
its strength being several times larger than FYH (TABLE VII) 
the new resultant force F’ determining the craft's roll 
behavior is directed well below the center of gravity 
(although F being directed above it). Thereby an outward 
roll angle is generated. 




Figure 3 - DEADRISE AND DRAG FORCES IN A PORT TURN 
(CG= CENTER OF GRAVITY, HL= WATER LINE) 
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Also from Fig. 3 the effect of considering the 
center of gravity versus the transom deadrise angle on the 
force generating the inward or outward roll moment can be 
seen. The center of gravity deadrise angle and its 
corresponding forces are shown in solid or dashed lines with 
F being directed above the center of gravity. The transom 
deadrise angle, shown in dotted lines, generates a force F 
being directed below the center of gravity producing an 
outward roll moment which gets reinforced by adding 
vectorially FYD. From this the general stability check ' 
yields a smaller outward roll moment and angle using the 
center of gravity versus the transom deadrise angle. This is 
in agreement with the results given in TABLE VI. 

Next, some simulation runs with arbitrary deadrise 
angles chosen such that small positive and negative steady 
state roll angles result are presented in order to, first, 
find out which deadrise angle would generate the moments 
necessary to simulate flat turns, and, second, verify the 
trend of the change in the respective moments which has been 
observed for the change from transom to center of gravity 
deadrise angle in TABLE VII. The value of the lateral added 
mass coefficient in these runs was 0.4 as it was in all runs 
listed in TABLE VII. Considering the roll moments over the 
range of deadrise angles from 78.8 to 28.7 degrees, the 
moments due to the bow seal, stern seal, propulsion and 
plenum pressure change consistently in magnitude while those 
due to sidewalls, rudder and aerodynamics reach some extreme 
value and then increase again. Regarding the steady state 
roll angles listed in TABLE VII, it can be seen that the 
arbitrary chosen deadrise angles happen to be nearly 
symmetric about the deadrise angle than would generate zero 
roll angle (flat turn) under these simulation conditions. 
Keeping this symmetry in mind one finds that the values of 
the roll moments due to bow seal (FKBS) , stern seal (FKSS) , 
propulsion (FKP) and plenum pressure (ABPB*PHI*Z) are also 
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symmetric with respect to zero moments. For propulsion and 
plenum pressure the sign of the moments are identical to 
that of the roll angle while for bow and stern seal moments 
the signs are reverse. Prom the results given in TABLE VII 
with positive steady state roll angles Figure 4 has been 
drawn where the dotted line represents the craft's initial 
state, the dashed line represents the craft's steady state 
considering the transom deadrise angle and the solid line 
represents the craft's steady state (smaller roll angle) 
considering the center of gravity deadrise angle for the 
vertical force calculation. Same markings apply to the 
moment vectors shown from which their change in magnitude 
for both deadrise angles can be seen. 




Figure 4 - ROLL MO.^ENTS IN A PORT TURN 
(ROLL MOMENTS AND ANGLE NOT TO SCALE) 



Legend: AED Aerodynamic 

ES Bowseal 
SS Sternseal 
RUD Rudder 
Deadrise angle: - 



CG Center of gravity 
PP Plenum pressure 
SW Sidewalls 
P Propulsion 
- - Transom 
Center of gravity 
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Sines the sinuilation run with the center of gravity 
deadrise angle being used for the vertical f-orce computation 
gave results that were mors favorable (better damped time 
histories and steady state values closer to those measured 
experimentally) this change in vertical force computation 
was implemented into the XR-3 Loads and Motion Program. 
Starting with the deadrise angle given in Table 7 for draft 
up to seven inches, again straight line approximations for 
the effective deadrise angle at the center of gravity 
location (station 14) for larger drafts have been performed 
resulting in an almost linear relationship. This part of 
subroutine SIDEWALL is now valid for any draft and contains 
the following statements which have been used in this form 
for all simulation runs to follow in this study: 

DDRAMG=0.0 

IF (DS.GT. 0.5833) DDRA NG= (DS-0 . 5 3 33) * 0 . 0629 
DRANG=1 .021+DDRANG-PM1 *PHI 
CTNDR=COTAN (DRANG) 

FZHOLD (J) =FZH (J) 

FZHDRP (J) =PM1 *FYH (J) *CT NDR* PROMO 1 
FZH (J) =FZH(J) +FZHDRP (J) 



# 

FK= (FZH (2) -FZH (1) ) * YS W + FKD-FY *Z S 

It should be noted that vertical projected forces 
from both sidewalls are used in computing the roll moments 
in the XR-3 Loads and Motions Program. Section III.E of this 
thesis will be concerned with the effect of changing the 
simulation program to consider the deadrise projected force 
only from the outward sidewall in a turn maneuver as it is 
discussed in Sect. III.B.2 . 
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c 



CROSS-FLOW DRAG 



In the preceding section the lateral force FYD whi 
due to cross-flow drag has been found to be rather 
compared to FIB, thereby changing the direction o 
resultant force significantly. FYD is comput 
subroutine SIDEWALL by suaning rhe contributions DF o 
stations : 



DF(I,J) = -HRH0*CDSW*VREL*A3S (VRSL) *DSWAV (I, J) 



where 



DF (I 


fJ) 


is the lateral drag 


contributed 


by 






the j-th element of 


each (i-th) 


sidewa 


HRHO 


is 


RHO/2. 






VREL 


is 


relative velocity 






DELX 


is 


the incremental dist 


ance used for the 




we 


tted draft calculati 


on along each side 


DSWA 


V is 


the average wetted 


draft over 


one 




incremental distance 






CDSW 


is 


the cross-flow drag 


coef f icient 


for on 




sidewall . 







The cross-flow drag coefficient being used in the p 
XR-3 Loads and Motions Program is the value correspond 
that for flew normal to a long flat plate being CDSW= 
Ref. 3 points out that an investigation of cross-flow 
forces and coefficients (before the vertical added ma 
had been reformulated as discussed in Sect. III.A.1) 
better agreement between theoretical and experi 
results if CDSW=2.f), representing the normal 
coefficient for a sharp flat plate, had been used in 
After the reformulation of A33 the agreement b 
modified theory and experiment had very much improv 



ch is 
large 
f the 
ed in 
f all 



wall 



resent 
ing to 
1. 28 . 

drag 
ss A 33 
showed 
mental 
drag 
stead, 
etween 
ed , so 
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that the former drag coefficient could be used again. Since 
in this thesis work rhe change in l33s did not effect the 
roll behavior significantly (TABLE IV ) , some turn maneuvers 
at 20 kn and 15 degrees rudder angle have been simulated for 
various drag coefficients in order to investigate which CDSW 
would better match simulated and experimental steady state 
roll angle. From TABLE VII a decrease in roll angle and 
roll moment is desireable and therefore the cross-flow drag 
force FYD should be decreased via reducing CDSW. The steady 
state values are shown in TABLE VIII with the lateral added 
mass coefficient being unchanged C =0.4. 



TABLE_VIII 

Stead y state co nditions for va rious CD SW 
2 0 kn t urn, 15 deg port rudde r angl e 



CDSW 


1 . 28 


1 . 16 


1 . 0 


0.7 


Pitch angle 


0.43 


0.43 


0 . 44 


0.46 


Roll angle 


1 . 55 


1 .48 


1 .38 


1.14 


Draft 


5.78 


5.80 


5.83 


5.91 


Speed 


19.39 


19.36 


19.31 


19.19 


FKBS 


-343.0 


-330.5 


-308. 1 


-252.8 


FKSS 


- 2.7 


- 2.6 


- 2.4 


- 2.0 


FKSW 


207.8 


175.0 


127.8 


40.0 


FKRUD 


- 82.2 


- 47.9 


- 2.3 


80.4 


FKAED 


- 56.6 


- 58.7 


- 61.6 


- 67.6 


FKP 


0 . 14 


0.14 


0. 13 


0.11 


AEPB*PHI*Z 


276.4 


264.6 


246 . 5 


201.9 


FYH (P,S) 


-4 3/-125 


-45/-125 


-48/-124 


-57/-1 21 


FYD 


-654.9 


-616.9 


-561 . 0 


-436.0 


Comparison 


, with Ref 


. 6 (Table 


VI) shows close agree 


in roll angle 


(1.360) 


for the run 


with CDSW=1 


.0 . But t 
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trial runs have been executed with arbitrary smaller drag 
coefficients because an increase in drag coefficient (e. g. 
CDSW=2.0 for sharp flat plate as in Ref. 3) would have 
enlarged the disagreement between simulated and experimental 
data. Also, from plot 29 for CDSW=1.28 and plot 31 for 
CDSW=1.16 it is seen that a smaller drag coefficient results 
in more damping in the roll response, in this case six 
versus nine cycles of distinguishable oscillations. Since 
the choice of CDSW=1.0 cannot be declared to be appropriate 
for the actual shape of the XR-3 sidewalls it will not be 
considered any further in this study. CDSW=1.16 generated a 
4 % change (about 15 % is desired) toward the experimentally 
measured roll angle and its effect will be investigated 
together with the present cross-flow drag coefficient 
CDSW=1.28 in the studies to follow. 

Shown below in Figure 5 is the dependence of steady 
state roll angle in turns generated by a 15 degrees rudder 
angle to port at 20 kn on the cross-flow drag coefficients 
used in the simulation runs (TABLE VIII) , 

Roll angle 




Figure 5 - S.3. ROLL ANGLE VERSUS CDSW 
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D. THRUST MAPPING 



In turn maneuvers as they have been discussed in the 
preceding sections another important point to be considered 
is the change in thrust on both engines. Since the engines 
mounted on the XR- 3 provide constant power for a given 
throttle setting, the actual thrust delivered during a turn 
maneuver may vary even though the throttle remains fixed. 
In all simulation runs shown in this thesis until now 

constant thrust has been used. In Ref. 6 experimentally 

measured reduction in thrust is given for turn maneuvers at 
20 kn and various rudder angles. From this and Ref. 7 
providing more detailed information the thrust-rudder 
dependence was as listed below. 

TABLE_IX 

R eduction in th rust for va ri ous rud der angles 
in turn maneuver s at 20 kn 



Rudder angle 
(deg) 



total reduction in thrust {Ref. 6) 

(%) 



5 

9 

12 



2.8 

3.3 

4.3 



TABLE IX has been used for thrust mapping (Block 16 of 
the data input) in the simulation runs. Due to 
unavailability of more precise data it has been assumed that 
the total reduction in thrust is contributed from both 
engines in equal amounts. But this may not be quite true, 
since due to the craft's roll angle a difference in 
immersion depth for both engines could result in unequal 
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i m 



I 



i 




changes in thrust. TABLE X shows the steady state values of 

these simulations considering two cross-flow drag 

coefficients, deadrise force contributions from both 

sidewalls computed at the center of gravity as well as both 

lateral added mass coefficients C =0.4 (0.8). 

h 



Steady 


state conditions usina 


thrust mapping 


in 


20 kn turns at 


various 


rudder angles 




Rudder angle 


Roll angle 


(deg) 


Speed 


(kn) 


(deg) 


simulation 


Ref. 6 


simulation 


Ref. 6 


CDSW=1 .28 


5 


0. 30(0.25) 


0.09 


19.3(19. 2) 


19. 3 


9 


0.71(0.61) 


0.28 


18.9(19.0) 


19.5 


12 


0.99(0.39) 


0.58 


18.5(18.6) 


19.2 


CDSW= 1 . 16 


5 


0. 28 (0. 23) 


0.09 


19.2 (19. 2) 


19.8 


9 


0.67(0.57) 


0.28 


18.9(19.0) 


19.5 


12 


0.94(0.34) 


0.58 


18.4 (18.6) 


19.2 



Comparing the steady state roll angles obtained if 
thrust mapping is used with those listed in TABLE VI (no 
thrust mapping, deadrise force computation at the center of 
gravity ) , a change in roll angle by 6 to 14 % for C =0.4 and 



by 14 to 20 


% for C =0.8 
h 


toward 


the measured 


angle 


can be 


observed . 


The steady 


state 


roll angles 


found 


in the 



simulation runs still differ from the roll angles given in 
Ref. 6 by a factor of about two and three but can be 
expected to be a little closer to these if proper thrust 
mapping, i.e. different maps for port and starboard engine, 
is used. By this also the steady state velocities which 
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differ by about 2.5 % can be expected to come in tetter 
agreement, i.e. to increase for smaller roll angles. 

Hegarding the roll angle time responses for simulation 
runs using both lateral added mass coefficients (0.4 and 
0.8) and 12 degrees rudder angle the following 

characteristics have been found: the runs using thrust 

mapping (Table X, plot 33(35)) show almost the same 

transient behavior as do the runs without thrust mapping 
(Table VI, plot 21 (23)) with about ten cycles of oscillation 
before the transients die out. But the runs using thrust 
mapping need quire a long time (t>45 sec) before they reach 
the steady state roll angle. 

The effect of using a larger lateral added mass 
coefficient (0.8) can be observed in plot 35. Comparing 
this with plot 33 , the roll angle response in plot 35 shows 
a smaller peak roll angle (1.25° versus 1 . 4 ? 0 ) and the 

transients have smaller amplitudes. But both die out after 
ten cycles. This effect has already been noted in Section 
III. B. 4. 



Investigating the effect of changing the coefficient 
CDSW=1.28 to 1.16 by comparison of the roll responses shown 
in plots 33 and 37 or 35 and 39, a significant damping 
effect as has been found previously in Section III.C can not 
be observed here. 
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E. 



DEADRISE FORCE OF OUTWARD SIDEWALL 



This part of the investigation on how to effect the roll 
behavior of the testcraft simulation is again concerned with 
the deadrise angle. 

In the present XR-3 Loads and Motions Program deadrise 
forces for both sidewalls are computed with the deadrise 
angle appropriately corrected by the roll angle for port and 
starboard sidewalls. The idea behind this, e.g. in a port 
turn, is to consider a pressure buildup (large FYH(S)) on 
the outboard side of the starboard sidewall and a reduction 
in pressure on the outboard side of the port sidewall (small 
FYH(P)). Thereby an upward vertical force FZHDRP (with 
negative sign) is obtained for starboard side, while a 
downward vertical force (positive sign) is computed for port 
side. Both vertical forces are then added to the buoyancy 
force of their side. The forces obtained by this approach 
are shown in Figure 6. 

In Ref. 4, however, a different philosophy in regarding 
the acting forces is presented. As discussed in 
Section III.E.2 an approximate check on the craft's roll 
stability can be made if the SES geometry and the vertical 
location of its center of gravity are known (Fig. 1) . In a 
port turn the principal force effecting the roll behavior 
arises from the outbo ard side of the starboard sidewall 
(relative high pressure on the structure there due to wave 
buildup) while there is only a small force (due to small 
pressure change) on the inboard side of the port sidewall 
contributing a roll moment only if there exists a roll 
angle. Now both vertical forces are directed upward. This 
approach and the resulting forces are shown in Figure 7. 
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Figure 6 - ACTING DEADHISE FORCES (TWO SIDEWALLS) 




Figure 7 - ACTING DEADRISE FORCES (ONE SIDEWALL) 



Note ; Since the vertical f 
sidewall is rather small, 
neglected in the roll stabil 



orce contributed by th 
as shown above, i 
ity check (Sect. III.B 



e inboard 
t can be 

. 2 ) . 
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To investigate the effect of this interpretation of the 
acting forces two statements controlling the computation of 
deadrise forces via the sign of the rudder angle RUDANG 
(port rudder is positive) and the sidewall under 

investigation (PM1=-1 for port side) have been added to the 
computation procedure for the vertical projected force ; 

DDEANG=0.0 

IF (DS.GT. 0.5833) DDRANG= (D S-0 . 58 33) *0 . 0629 
DRANG=1 .021 +D DRANG- PM1* PHI 
CTNDR= COTAN (DRANG) 

ROD SIG=SIGN ( 1 . , RUDANG) 

IF (RODSIG.NE. PM1) CTND R= PHI ^i'T AN (PHI) 

FZHOLD (J) =FZH (J) 

FZHDRP (J) =PM1 *FYH (J) *CT NDR* PROMO 1 
FZH (J) =FZH (J) +FZHDRP (J) 

The steady state values for simulated turn maneuvers using 
thrust mapping and center of gravity deadrise angle are 
shown below considering two cross-flow drag coefficients and 
both lateral added mass coefficients C =0.4 (0.8). 



TA BLE XI 



Steady state conditions 



in_2^ 


kn turns at 


various 


rudder angles 




Rudder angle 


Roll angle 


(deg) 


Speed 


(kn) 


(deg) 


simulation 


Ref . 6 


simulation 


Ref. 


CDSW=1 . 28 


5 


0.43 (0. 52) 


0.09 


19.4 (19. 5) 


19.3 


9 


0.90(1 .00) 


0.28 


19. 1 (19. 4) 


19.5 


12 


1- 22( 1. 34) 


0.61 


18.7 (19.0) 


19. 2 


CDSW=1 .16 


5 


0.4 1(0.5 1) 


0.09 


19.3(19.5) 


19.3 


9 


0.87(0.98) 


0.28 


19.1 (19. 3) 


19.5 


12 


1.17(1.30) 


0.61 


18.6(19.0) 


19.2 



6 
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Comparing the steady state values for simulated 
tables X and XI , the effect of using only th 
sidewall for the deadrise force computation can be 
to give less favorable (little larger) steady s 
angles but mere favorable (larger) steady state ve 
This is exactly what has been expected since, as 
Figures 6 and 7, if deadrise force contributions f 
sidewalls are used, the deadrise force from 
sidewall will reinforce the one from the starboard 
toward a less outward roll moment (Fig. 6) while 
the philosophy depicted in Fig. 7 the vertical fo 
the port sidewall counteracts the starboard deadri 

Regarding the roil angle responses shown in ip 

for 12 degrees rudder angle with CD3W=1.28 (plots 4 

and CDSW=1.16 (plots 45 and 47) and comparing these 

corresponding ones from the previous section (pi 

and 37, 39) the change in the vertical force co 

considering the outward sidewall only is seen to 

about 8 versus 10 cycles of transient, oscillations 

amplitude. Using C =0.8 (plots 43 and 47) ha 

h 

significant damping effect and makes the roll angle 



more resemble a step, which is in accordance 
experience gained by Ref. 7. Also, comparing the 
with plots 35 and 39, the negative roll angle occ 
the moment when the rudder has been introduced does 
up any more. 
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VBHTICAL LOCATION OF CENTER OF GRAVITY 



In Section III.B.2 it wa 
testcraft's roll bahavior 
geometry and the location o 
known. This section will 
different vertical locations 
simulated XR-3 roll behavior 



s discussed how to check for the 
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turn provided the 
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center of gravity 


are 
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cone 


;erned with the effect 


of 


of 


the 


center of gravity on 


the 



From experimental measurement the location of t 
of gravity of the XR-3 has been determined to be 
forward of the transom and Leo and Boncal establi 
vertical location at 2.54 ft above the keel 
longitudinal center line of the craft. As modifica 
the testcraft were introduced (e.g. engines and s 
horizontal location of the center of gravity w 
established experimentally. No such changes were ma 
vertical location. In order to investigate the se 
of the steady state roll angle to vertical locatio 
center of gravity, reductions in height (in accorda 
the stability check geometry depicted in Figu 
increments of 0. 1 ft were entered into the Loads a 
Program. 
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Simulation runs with constant thrust and initial speed 
of 20 kn introducing a 15 degrees port rudder angle, 
deadrise force contributions from both sidewalls computed at 
the center of ^ gravity, the sidewall drag coefficient 
CDSW=1.28 and the lateral added mass coefficient being 0.4 
were executed and the results are shown in Tanle XII. From 
the listed data it follows that with a lower vertical 
location of the center of gravity a tendency toward an 
inward roll angle (here actually a smaller outward roll 
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angle) can be achieved. As shown, a decrease in vertical 
distance by 0.1 ft resulted in about 0.1 degree decrease for 
both peak and steady state roll angle (almost linearly 
related) . From the corresponding plots as indicated in 
Table XII a roll damping effect due to the relocation of the 
center of gravity can hardly be recognized. The number of 
cycles of transient oscillations and the steady state speed 
were about the same for ail runs listed in Table XII. 



TABLE_XII 

Roll conditio ns in 2 0 kn r ur n with con stant thru st. 
a nd 15 degr ee s p o rt r udder angle 



vertical location 


of CG 


experimental 


2.54 


2.44 


2.34 


2.24 


testcraft data 


1.55 


1.43 


1.32 


1 .22 


1.36 


1.67 


1.53 


1 . 46 


1 .35 




2. 13 


2. 07 


1 .98 


1.91 




8 


8 


8 


8 




13-14 


49-50 


51-52 


53-54 




19.4 


19.3 


19.3 


19.3 


18. 7 


shows 


the 


effect 


of ' 


changing the vertical 



steady state 
roll angle 
peak 

roll angle 
peak 

roll rate 
cycles of 

transients 
refer to 
plots 

steady state 
speed 



location of the center of gravity for simulation runs at 
20 kn with 12 degrees port rudder angle and thrust mapping 
applied, deadrise force computation for both sidewalls at 
the center of gravity and the coefficients being 0.4 for the 
lateral added mass and CDSW=1.28 for the sidewall drag. 
These conditions were similar to those used in Table X. 
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From Table XIII the effect of reducing the vertical location 
of the center of gravity in 0.1 ft increments over a 
reasonable range can be observed to also result in linear 
reductions in maximum and steady state roll angle for the 
case if thrust mapping is used. Nonlinear reductions appear 
for the peak roll rate. The number of cycles during the 
transient period (ten) is again not effected by the vertical 
relocation of the center of gravity. Regarding the 
corresponding plots a very slight damping effect can be 
observed for reduced vertical locations. The results given 
in Table XIII are graphically represented in Figure 8 . 



TA BLE XIII 

Roll con dit ions i n 20 _k n _ t u r n _w i t h _ t t ma pp in g 
a nd 1 2 degre e s port rudder angle 



vertical location of CG experimental 

2.54 2.44 2.34 2.24 testcraft data 



steady state 

roll angle 0.99 0.91 0.34 0.77 



0. 53 



pea k 

roll angle 1.47 1.36 1.25 1.16 



peak 

roll rate 2.14 1.92 1.74 1.58 



cycles of 

transients 10 10 10 10 



refer to 
plots 

steady state 
speed 



33-34 55-56 57-58 59-60 



18.5 18.4 18.4 18.4 



19.2 
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Figure 8 - ROLL CONDITIONS VERSUS 
VERTICAL LOCATION OF CENTER OF GRAVITY (TABLE XIII) 
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ROLL EAMFING DOE TO VERTICAL WAVE GENERATION 



One part of subroutine SIDEWALL is co 
roll damping effect which is due to vertical 
during a roll motion. The computation of 
term was developed by Oceanics Inc. and 
program by Menzel [Ref. 2]. The princi 
development cf this addition were not provid 
Inc./ however, the apparent idea is to red 
roll moment FK by some value which is de 
craft's roll rate and its draft, thereby 
added mass A33. During the study of the e 
additional terms a discrepancy was noted i 
used in the calculations. 
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The roll moment FK is calculated by' the equation given 
in SDWL1950 and considering the expressions leading to this 
statement it is found that 



dim [FK] = lb . ft 

which is the correcx dimension for a moment vector. Now 
considering the dimension of the roll damping terra 

dim [ PR0M02^YSW*YSW*BC2*P/PI ] 

which is subtracted from FK (SDWL2200) it is found 



PRCB02 and 


PI are 


dimensionles 


dim 


[ 


YSW ] = 


ft 








dim 


C 


P] = 


rad/s 


ec 




dim 


[ 


BC2] = 


dim 


L 


AC2] 










dim 


c 


FC2 * 


length ] 






= 


dim 


c 


F2 * 


cos(x) * 






= 


dim 


c 


A33 * 


length ] 



length ] 
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= dim [ RHO B * B ♦ length] 

= (lb . sec2 /ft*) . ft . ft . ft . 

So for the roll damping term 
dim [YSW * ysw * P * BC2] 

= ft2 . (rad/sec) . (lb . sec2/ft*) . ft^ 

= lb . ft . sec 

which is not the correct dimension for a moment expression. 
Assuming that only the given terras may be involved in the 
roll damping calculation, it is possible to come up with the 

proper dimension if the roll rate P were squared but keeping 

its sign, thus getting 

roll damping term = PR0i^02*YSW*YS» *BC2*P*ABS (P) /PI . 

Also the vertical added mass should be associated with some 
‘velocity squared* - expression which is absent in the 
original version of this part of the program but could be 
generated by squaring the roll rate P. The effect of the 
modified expression will be an increase in the damping 
effect for roll rates P>1.0 which exist only in the initial 
phase of the turn maneuver (maximum roll rate in all 
simulation runs was about 2.0) . But for most of the run 

length P is less than unity and the supposed modification 
will decrease the damping effect especially when approaching 
steady srate (P=0.0) . 

The PS0M02 - term provides a means to arbitrarily set a 
roll damping factor which, as stated in Reference 2, has 
been found from experiments to be 16.0 . Until now 

PROMO 2=1.0 has been used in this thesis work. So the effect 
of changing it to the experimental value is investigated 
next. In the runs to be discussed turn maneuvers at 20 kn 
with a 12 degrees port rudder angle have been simulated with 
deadrise forces being computed at the center of gravity, 
using the experimental damping factor PROMO2=16.0 and the 
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coefficients CDSW=1.28 for the sidewall drag and C =0.4 for 

h 

the lateral added mass. 



1. Thrust is held constant with deadrise force cont 
from both sidewalls 

Comparing the obtained plots 61-62 with plot 
(PfiOMO2=1.0) only a slight damping effect on r 
response can be observed (peak roll angle 1.45 
1.47°), but the number of transient oscillations 
from eleven to nine. The peak roll rate decreased f 
to 2.10 for PROMO2=16.0 and the damping effec 
pronounced in the roll rate time history. 



ributions 

s 21-22 
oil angle 
0 versus 
decreases 
rom 2.14 
t is more 



2. Thrust is mapped with deadrise force contributions from 
both sidewalls 

Comparing the plots 63-64 with the corresponding ones 
(33-34) for PROM02=1.0 , again only a slight damping effect 
can be noticed. The peak roll angle is again reduced from 
1.47 to 1.450 and the number of cycles of transient 
oscillations is uneffected. 



If the damping term is changed to include the P*A3S(P) 
expression (plots 65-66) there result eleven versus nine 
cycles of transient oscillations with little higher 
amplitudes (peak roll angle is again 1.47o as it was for 
PROBO2=1.0) . The roll rate time history also shows little 
larger amplitudes (peak roll rate is 2.14 versus 2.10) . 



3. Thrust is mapped with deadrise force contribution from 

the outward sidewall only and lateral added mass 

coefficient C =0.8 
h 

Comparing these plots (67-68) for PROMO2=16.0 with plots 
43-44 where PRCMO2=1.0 again only a slight damping effect 
in roll angle response can be noticed with the number of 
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transient oscillations being uneffected. Using the P*ABS(P) 
term the slight damping effect shown in plot 67 is cancelled 
as shown in plot 69. 

The above observations show that the experimentally 
determined roll damping factor PROMO2=16.0 has only a 
negligible damping effect. Considering the roll damping term 
again as it was added to the present XR-3 Loads and Motions 
Program in 1974, there was a 1/PI - factor which actually 
cancelled the PI - factor contained in the former expression 
for the vertical added mass A33S (Section III. A). Since the 
PI - factor already has been eliminated by correcting the 
A33S - expression, there is no further^need to cancel it in 
the damping term by the 1/PI - term which therefore also 
could be eliminated. This change is expected to cause an 
increase in sidewall roll damping by a factor of PI and has 
been investigated by an additional run. 

4. Thrust is mapped with deadrise force contribution from 
the outward sidewall only, roll damping term being 
16 .0*YSW*YSW*BC2*P 

and lateral added mass coefficient C =0.8 

h 

From plots 71-72 again only a slight damping effect versus 
plots 67-68 can be noticed. 

The investigation described in this section shows that 
the roll damping due to vertical wave generation as included 
in subroutine SIDEWALL is not very effective even if the 
factor PR0M02 in the damping term is increased from 16.0 to 
50.0 . There also is a laclc in matching dimensions 
(statement SDWL2200) . 
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IV. PROPULSION AND RUDDER SUBROUTINES 



In all simulation runs executed during this study the 
XR-3 Loads and Motions Program as given in Ref. 2 has been 
used except for the modifications mentioned in this work. 
This simulation program contains Forbes* version [Ref. 9] 
of 'SUBROUTINE PROP' which included some revisions of the 
original version given by Leo and Boncal [Ref. 8]. Since in 
TABLE VII (Section III.B.U) the roll moment contributed by 
the propulsion system (FKP) is rather small (0.1 to 0.2) an 
additional run using Leo and Boncal' s version has been 
executed under the conditions as stated for TABLE VII 
(center of gravity deadrise angle). The steady state values 
are shown below (in parentheses are the values obtained 
using Forbes' version) : 



pitch angle 


0. 44 (0 .45) 


draft 


5. 73 ( 


5.81) 


roll angle 


1. 67 (1 .53) 


speed 


19.39 (19.37) 


FYSW 


-909. 4 (-821 . 9) 


FYRUD 


- 54.7 ( 


32.7) 


FYAED 


- 2 1 . 3 (- 20.0) 


R*V*AM 


888.4 ( 


809.2) 


FYP 


97. 0( 0.0) 








FKBS 


-348.0 (-336 .4) 


FKSS 


'- 3 . 0 (- 


2.7) 


FKSW 


269. 9 ( 2 12. 1) 


FKRUD 


150.3 (- 


89.9) 


FKAED 


- 63. 5 (- 56 .0) 


FKP 


-305.2 ( 


0.14) 


A3PB*PHI*Z 


299. 5 ( 273.5) 








FYD 


-729. 2(-651. 4) 


FYH (P/S) -43/ 


-137 (- 
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Regarding these values drastic differences are found for 
FYP and FYRUD as well as for FKP and FKEUD, though the 
steady state values for pitch and roll angle, draft and 
speed are close to those obtained using Forbes* version. 
The pitch and roll angle responses are given in plots 25-28. 
Comparing these with the corresponding ones using Forbes’ 
version (plots 13-16) it is seen that Leo and Boncal's 
version generates more oscillatory roll and pitch angle 
responses. Therefore a review of subroutine PROP especially 
the computation of the respective forces and moments has 
been carried out and are discussed below. 

The review resulted in force and moment equations which 

are identical to those given by Leo and Boncal [Ref. 8]. 

Another point that supports the supposition to use this 

\ 

version of subroutine PROP is the fact that it provides a 
reasonable roll moment due to rhe propulsion system (Forbes' 
version provides almost zero roll moment) . In order to test 
the effect of changing subroutine PROP to its original form, 
simulation runs 1. and 2. from the preceding section have 
been repeated and their time histories for the roll behavior 
are shown in plots 73-74 and 75-76, respectively. Comparing 
these with the corresponding ones using Forbes* version of 
subroutine PROP, plots 61-62 and 63-64, it is observed that 
the sready state and peak roll angles are larger by about 
7 % and the number of transient oscillations is uneffected 
by this change. 

Next a review of subroutine RUDDER has been carried out 
in order to investigate the cause of the change in sign for 
FYRUD and FKRUD depending on the propulsion subroutine used 
(Forbes* or Leo and Boncal's version) although identical 
RUDDER-subroutines have been used in both cases. The 
important equations to be considered are : 
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DSR=Z+ZS-XS*THETA 
ENDFAC=1. + DSR/ (DSR + RSPAN) 

VH=V + XR*R - ZR*P 
EFFANG=RUDANG - ENDFAC*VH/CJ 
FY=RHO*U*U*RAREA*ENDFAC*RCLB*EFFANG 
• 

FK= -ZR*FY 

The above equations compute 

- draft at the rudder location 

- endfactor depending on the rudder length below the 
craft's keel 

- sideslip velocity 

- effective rudder angle depending on the craft's forward 
and sideslip velocities 

- lateral force depending on rudder area and effective 
rudder angle 

- roll moment due to rudder . 

The sign of the roll moment FK is dependent on the sign 
of the lateral force FY which depends on the sign of the 
effective rudder angle EFFANG. This in turn depends of 
course on the introduced rudder angle, but the magnitude of 
the second term of this equation (sideslip angle) is 
responsible for the actual sign of EFFANG. VH/D represents 
the tangent of the sideslip angle which - applying small 
angle approximation - is subtracted from the introduced 
rudder angle (see Figure 9). Since ENDFAC=1.73 and 0 was 
the same for both runs, a difference in VH which was larger 
in Leo and Boncal's version by about 50 % caused the 
negative sign (versus a positive sign using Forbes' version) 
for the effective rudder angle EFFANG and thereby FY. So 
the only way to introduce a change toward an inward roll 
behavior here (FKRUD was positive) is to first check the 
input value for RSPAN and second to eventually modify the 
craft's engines to have larger immersion into the water. 
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with regard to the first item, if RSPAN would increase, 
ENDFAC would decrease as would the magnitudes of SFFAHG 
(negative), FY (negative) and FK (positive). 




RUDANG rudder angle 

EFFANG effective rudder angle 

Figure 9 - EFFECTIVE RUDDER ANGLE IN A TORN MANEUVER 



Considering next the fact that 
had been replaced in 1976 a recent 
and propeller location yielded the 
input, values shown in parentheses) 



the engines of 
measurement of 
following data 



the XR-3 
the rudder 
(or igina 1 



XPO = -1 . 50 (- 1. 275) 

ZPO = -1 . 333 (-0.604) 

XRO = -1 . 083 (- 1. 125) 
RSPAN= 1.333 ( 1.21 ) 
RARSA= 1.42 ( 0.68 ) 
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Two of the runs discussed in the preceding section have 
been repeated using the original roll damping term and the 
new propulsion and rudder input parameters. 

1. Thrust is mapped with deadrise force contributions from 

both sidewalls with C =0.4 . 

h 

Comparing plots 77-78 with the corresponding ones for the 
former rudder and propulsion parameters (plots 63-64) the 
expected effect due to a decrease in FKRUD can be observed 
resulting in smaller peak roll angle (1.34° versus 1.450) 
and a significantly decreased steady state roll angle (0.65o 
versus 0.99® at t=50 sec) . The transient period of the roll 
angle response in plot 77 shows the same number of transient 
oscillations as plot 63 but the response did not reach 
steady state even after 45 sec . The steady state roll 
angle could be estimated to be about 0.3 to 0.4 degrees. 

2. Same conditions as in the previous run but with deadrise 

force contribution from the outward sidewall only and 

C =0.8 . 
h 

Regarding plots 79-80 versus plots 71-72 a similar effect as 
mentioned before is observed. Now the roll angle response 
shows a more oscillatory transient period. The peak roll 
angle reduced from 1.46 to 1.36 degrees while the roll angle 
at t=50 sec changed from 1.34 to 0.97 degrees. For this run 
the steady state roll angle could be estimated to be about 
0.8 degrees since the roll angle response again did not 
reach steady state for the run using the recently determined 
rudder and propulsion parameters. The initial small 
negative roll angle which appeared in the previous run 
(plot 77) at the moment when the rudder morion is 
introduced, does not show up any more. The overall roll 
response shown in plot 79 does not fall off as steeply 
(compared to plot 77) if, in a turn maneuver, deadrise force 
contribution from the outward sidewall only is considered. 
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From these runs it is obvious that the rudder and 
propulsion parameters of the new engines result in a more 
favorable XR-3 roll behavior with a significantly less 
steady state outward roll angle. At the time of this 
sensitivity study Reference 7 pointed out that there 
actually is no roll angle indicated by the measuring devices 
in a port turn at 20 knots and 12 or 15 degrees rudder 
angle. Keeping the device sensitivity of ±0.5 degrees in 
mind, it can be concluded that the steady state roll angle 
provided by the present XR-3 Loads and Motions Program 
including the before mentioned program modifications by the 
author, is quite well in agreement with the actual craft 
data. Further efforts should be undertaken to obtain a 
better damping effect during the transient period. 

Remark : 

Although the new engines extend about 0.5 ft deeper into the 
water than the old engines, the new input value for RSPAN is 
close to the former value. In Reference 1 the denominator of 
the second term in the equation for ENDFAC is defined to be 
the distance from the water surface to the bottom tip of the 
rudder. The equation for ENDFAC used in the present program 
accounts for changes in draft DSR in both the numerator and 
the denominator. The RSPAN-term should then be defined to 
be the distance from the craft's keel to the bottom tip of 
the rudder. The author suspects that formerly there has been 
chosen too large a value for RSPAN in the case of the old 
engines (RSPAN=1.21 ft, probably from water surface to 
bottom tip of the rudder) . Since the bottom tip of the 
rudder and the center of the propeller were located at 
almost the same distance from the craft's keel the 
magnitudes of ZPO and RSPAN should be about the same which 
is not true for the former chosen values. From these 
observations the author concludes that until now a too 
optimistic roll behavior resulted due to a too large value 
for RSPAN. 
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V. 



CONCLUSIONS AND R2C0MMENDATI0NS 



The preceding sections of this chapter reflect the 
results of an investigation on how the XR-3 roll behavior in 
turn maneuvers at 20 kn simulated with the Loads and Motions 
Program is affected by certain changes in parameters and 
forces. The results of this sensitivity study are 
summarized as follows : 

* Added mass effect 

A reformulation of the vertical added mass computation which 
reduced the A33s-value by a factor of 1/PI had only a little 
effect in changing steady state draft, pitch angle and 
required thrust for constant speed in straight runs. In 
turn maneuvers this change resulted in less damped responses 
for pitch and roll angle and affected the steady state 
values for pitch angle by +8.3 % , roll angle by -1.0 % and 
draft by +1.4 % . Considering the roll behavior, the effect 
of reducing the A33s-value by a factor of 1/PI in the 
vertical added mass computation was negligible. 

A change in the lateral added mass coefficient C =0.4 to 

h 

0.8 had the effect of reducing the amplitudes of the 

transient oscillations and the steady state roll angle by 

0.1 degree toward the experimentally measured value. Since 

in a turn maneuver with 15 degrees rudder angle an unstable 

craft behavior showed up which had not been experienced in 

practice, the value of C =0.8 is not realistic if deadrise 

h 

force contributions from both sidewalls are considered. If 
deadrise force contribution from the outward sidewall only 
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is used in the simulation 
added mass coefficient cause 
state roll angle by abou 
roll angle response is q 
resembles a step as it doe 
larger coefficient (0.8) see 
case of the deadrise for 
sidewall only. 

* Deadrise angle effect 

The deadrise angles over th 
sidewalls are not uniform. C 
is used in the deadrise fore 
(78.80) to the center of 
damped responses for pitc 
maneuvers. For all chosen 
turns at 20 kn the steady st 
by 23 to 38 % toward the c 
experimentally and documente 



program the increase in lateral 
s an increase in the steady 
t 0. 1 degree. But since now the 
uite significantly damped and 
s in practice, the choice of the 
ms to be appropriate for the 
ce contributed from the outward 



e length of the curved XR-3 
hanging the deadrise angle which 
e computation from the transom 
gravity (58.5°) resulted in more 
h and roll motion in turn 
rudder angles generating port 
ate roll angle changed favorably 
orresponding roll angle measured 
d in Ref. 6. 



* Cross-flow drag coefficient 

Changing the cross-flow drag coefficient (CDSW=1.28 for a 
long flat plate) to lower values, e.g. CDSW=1.00 , the 
testcraft's simulated roll behavior could be forced to match 
the measured steady state values. But since this 
coefficient, has been chosen arbitrarily and cannot be 
declared to be appropriate for the actual XR-3 sidewalls 
(CDSW=1.0 means that there is no effective drag coefficient) 
this CDSW-value has not been considered any longer. Another 
suggested cross-flow drag coefficient CDSW=1.16 [Ref. 3] 
resulted in about 5 % reduction in steady state roll angle 
and a shorter transient period (about 20 %) with smaller 
amplitudes of oscillation than did tjje runs using CDSW = 1.28. 
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* Thrust mapping 



i 

Performing thrust mapping in the simulation runs using the 
data provided in Ref. 6 and assuming that both engines lose 
equal amounts of thrust in turn maneuvers, the simulated 
steady state roll angles could be shown to agree better with 
the measured angles by about another 10 % but the roll angle 
response did not quite reach steady state even after 45 
seconds . 



* Deadrise force from outward sidewall only 

The simulation runs executed with the deadrise force 

computed for the outward sidewall only and thrust mapping as 

done before, showed less agreement for both C =0.4 (0.3) 

h 

with the corresponding measured steady state roll angles 

than did the runs considering deadrise forces from both 

sidewalls. But for C =0.3 the roll angle response was 

h 

significantly damped and looked more like a step, as was 
experienced in practice. 

* Vertical location of the center of gravity 

The given value for the vertical location of the center of 
gravity (2.54 ft above the keel) might not be quite correct 
due to the difficulties in experimentally determining it. 
The significant effect of the vertical center of gravity 
location on the XR-3 roll behavior is obvious from Figures 1 
and 8 and the discussion in Section III.?. Since it controls 
exceedingly the XR-3 roll behavior it is important that 
future investigators obtain an accurate measurement of the 
C.G. location before taking testcraft data. 



* Subroutine PROP 



A review of the propulsion subroutine revealed tha 
force and moment calculations as given by Leo and 
(Ref. 8) should be used instead of the version stat 
Forbes (Ref. 9). This change results in a 7 % larger o 
steady state roll angle but it provides a reasonable 
moment due to the propulsion system. 

* Roll damping due to vertical wave generation 

The roll damping calculation included in subroutine SI 
has only a negligible effect on the craft's net 
behavior over the range of the damping 
1 . 0<PROH02<50 . 0 . Since there also is a discrepancy i 
dimension of the roll damping term the source of thi 
of the program should be investigated. 

* Suggested areas for further investigation 



During the studies undertaken with the XR-3 Loads and 
Program the author found the following points either 
input data from the testcraft or in the program that r 
further investigation : 



■ Rudder/thr ust angle 

As learned from Ref. 7 the rudder/thrust motion is di 
introduced to the port engine while the starboard 
follows this via a metal bar connecting both engine 
starboard rudder and thrust vector may be off by 2 
degrees to either side, depending on the direction 
rudder/t nrust action. Thereby the roll behavior 
certainly be effected. Since rudder/thrust mappi 
available in the Loads and Motion Program, it is imp 
that accurate data be collected for each engine du 
test run. 
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■ XE-3 weight 

During each test run there is a possibility that the weight 
is not exact since, as learned from Ref. 7, the testcraft 
takes on water in the seals and thereby its weight increases 
causing larger draft which also effects the roll behavior. 
This additional weight can be substantial and attempts 
should be made to measure or estimate the additional weight 
to be used in the Loads and Motion Program. 

• Cross- flow drag coefficient 

Further efforts should be undertaken to establish the proper 
value of the cross-flow drag coefficient since CDSW=1.28 
might not be quite appropriate for the actual shape of the 
XR-3 sidewalls. 

■ Sideslip velocity 

Those parts of the Loads and Motions Program where terms 
contributing to the lateral velocity V are computed should 
be reviewed since a large value of V may result in a change 
in the effective rudder angle from positive to negative. 

In this sensitivity study it has been shown that the 
most significant effect on the simulated XR-3 roll behavior 
are contributed by 

* deadrise angle 

* vertical location of the center of gravity 

* rudder and propeller location . 

Less significant effects are contributed by 

* thrust mapping 



* cross-flow drag coefficient . 




k 



For further studies concerned with the XR-3 Loads and 
Motions Program for improved roll behavior representation it 
is recommended to use 



- the reformulated equation for the vertical added mass 
A33 

- deadrise angle at the center of gravity location 

- deadrise force contribution from the outward sidewall 
only 

- subroutine PROP as given by Leo and Boncal and in 
Appendix C of 'this thesis 

- thrust and rudder mapping. 
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APPENDIX A 



PLOTS 



In the table given below are listed all plots with the 
respective parameters used in the corresponding simulation 
run. The table contains the following abbreviations : 
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for applied parameters see first page of tnis appendix 
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for applied parameters see first page of this appendix 
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APPENDIX B 



MODIFICATIONS OF THE SIMULATION PROGRAM 



_ During the course of studies undertaken with the XR-3 
Loads and Motion Program some statements have been added or 
changed in order to improve the work with the simulation 
program, 

* MAIN program 
The statements 



have been modified and added to the MAIN-program, The 



purpose of the second statement is to start the writing of 
the output values desired if computer time has reached the 
demanded time (jobcard) reduced by IQUIT which is the 
expected amount of time (in 10 ♦ seconds) required to print 
the output. The user may input the desired value of IQUIT as 
the fourth parameter (110-Format, columns 36-45) on the 107 
data card. Default value for IQUIT is 600000 (60 seconds), 

* Subroutine INCON 

The added statements are 



COMMON /AIR/ PINF,RHOINF, GAM, IQUIT 
IF (MYTIME(DUM) ,LT. IQUIT) GO TO 12 



(MAIN0050) 
(HAIN066 1) 



IQUIT=600000 



(INCN0831) 



to set the default value and 



READ (5,204 1) IQUIT 



(INCN1455) 
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to read the set value for IQUIT thereby overriding the 
default value. 

* Subroutine RHS 

In Wenzel’s version [Ref. 2] of the XR-3 Loads and Motions 
Program the entered craft’s velocity 00 (kn) is changed in 
Subroutine INCON to 0 (ft/sec) by 

0 = 00*1.6889 (INCN4580) 

and in Subroutine RHS transformed back to units in knots 
(RHS 1850) by 

VEL = 0.5925*0 . (RHS 1850) 

From these two equations follows 
VEL = 0.5925*00*1.6889 

and for 00=20 kn there results VEL = 20.0135 kn . But since 
it is desired that VEL=U0 the statement has been corrected 
to 

VEL = 0/1 .6889 . 

Statements MAIN0840, RHS 1140 have been changed accordingly. 

* Subroutine INTGRL 

To this subroutine the statement 

CALL COLFIL (INT 1021) 

has been added in order to provide the output values already 
calculated to the user if the minimum time interval allowed 
in the integration process is undergone. Before this 
addition only the warning message 

’DELTA TIME LESS THAN 1.0E-6 - - JOB STOPS' 

appeared and the output values have been lost. 
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XR-3 LOADS AND MOTIONS PROGRAM 
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AND 2) INITIALIZE VARIABLES FOR CALCS 
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